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Abstract 



The Einstein-Boltzmann system is studied, with particular attention to the non-negativity 
of the solution of the Boltzmann equation. A new parametrization of post-collisional momenta 

fj ' , in general relativity is introduced and then used to simplify the conditions on the collision 

cross-section given by Bancel and Choquet-Bruhat in 3 . The non-negativity of solutions of 
the Boltzmann equation on a given curved spacetime has been studied by Bichteler [3] and 

hf\. Tadmon [20] . By examining to what extent the results of these authors apply in the framework 

of Bancel and Choquet-Bruhat, the non-negativity problem for the Einstein-Boltzmann system 
is resolved for a certain class of scattering kernels. It is emphasized that it is a challenge to 
extend the existing theory of the Cauchy problem for the Einstein-Boltzmann system so as to 
include scattering kernels which are physically well-motivated. 
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1 Introduction 

To understand how the gravitational field and matter interact in general relativity it is necessary to 
specify a suitable matter model. One frequently used type of model comes from kinetic theory, where 
matter is considered as a collection of particles described statistically. The Vlasov and Boltzmann 
equations are the equations describing this type of model and coupling to the gravitational field gives 
the Einstein- Vlasov and the Einstein-Boltzmann systems respectively. We refer to [T] and [TB] and 
their references for basic information about the Einstein- Vlasov system and the Einstein equations 
coupled to many different matter models. The present paper considers only the Einstein-Boltzmann 



system. Compared to the Einstein- Vlasov system, the Einstein-Boltzmann (EB) system has not been 
studied much in the decades since local existence was proven by Bancel and Choquet-Bruhat [21 [3] . 
Recently several further results have been obtained for the EB system. Mucha [TTJ [TJ] showed that 
the initial regularity conditions given in 2,3 can be weakened. Noutchegueme and Dongo [T3] and 
Noutchegueme and Takou [15] proved some global existence results for certain classes of spacetimes 
with symmetries. The relativistic Boltzmann equation on a given curved spacetime was studied in 
[HOED]- As for the nonrelativistic or special relativistic Boltzmann equations, one can find plenty 
of references, and we only refer to [5j El |8] . 

The EB system is a system of equations where the Einstein equations and the Boltzmann equation 
are coupled to each other. In harmonic coordinates, the Einstein equations take the form of the 
system of quasilinear wave equations 

- -9 lS f 9 1 i + H af} = T afi - -(g-< s T^)g a P, (1.1) 

for a Lorentzian metric g a p, where H a P is a rational function of g a p and d 1 g a p with denominator 
a power of the determinant \g\, and T a ^ is the stress-energy tensor of matter described by the 
Boltzmann equation. The Boltzmann equation describes the dynamics of a relativistic gas, the 
particles of which interact through binary collisions. The unknown is the distribution function / 
defined on the phase space, which is the tangent bundle T(M) of a four-dimensional manifold M with 
a metric g a p. If we assume that all the particles considered have the same mass, then momentum 
is confined to a submanifold P(M) C T(M), which is defined by g a pp a p = — 1, when the mass 
is normalized to unity. This submanifold is called the mass hyperboloid. As a consequence of this 
assumption the distribution function turns out to be a function of seven variables. The Boltzmann 
equation takes the form 

pa ^- r > ap0 w = QifJ) > (L2) 

where T l a o are the Christoffcl symbols of the metric g a p, and the collision operator Q, 

Q(f,f)= f f S(x,p,q,n)(f(p')f(q')-f(p)f(qj)dn^dq, (1.3) 

will be studied in detail in the next section. Here, for simplicity we abbreviate f(x,p') as f(p'), 
f(x,q') as /(<?'), and so on. The stress-energy tensor T a/3 is defined as 

«0 - f fl^ ^«<wM 



T ap = / /(a:,p)pV— <fc, (1.4) 

Jr3 -p 

where the integration over M 3 is understood as integration on the mass hyperboloid at x. The 
equations (|TTTj) - ([P)) form the EB system [71 IT7] . 

In this paper, we are interested in the strong solutions of the EB system constructed in [5] 
by Bancel and Choquet-Bruhat. In Section 2, the Boltzmann equation in a curved spacetime is 
studied and an alternative formulation using an orthonormal frame is explained. (Its relevance is 
discussed in Section 4.) Returning to the equations written in a coordinate frame as in [3] we 
introduce a new parametrization of the post-collisional momenta. Unlike in the nonrelativistic case, 
the post-collisional momenta can usefully be parametrized in many different ways in the relativistic 
case HH1 [19]. The parametrization in Section 2.1 can be thought of as a generalization of that 
in [S]. With this parametrization, we consider the original result of Bancel and Choquet-Bruhat, 
and simplify some conditions in their paper. The distribution function in the Boltzmann equation 
should be non-negative because of its interpretation as a number density. In mathematical terms 
this means that a solution / of the Boltzmann equation arising from non-negative initial data 



should itself be non- negative. This important property was not proved in |3|. A central aim of 
the present paper is to clarify this question, taking into the account the available results in the 
literature. This non-negativity problem for the EB system is considered in Section 3. In the case of 
the Boltzmann equation on a given curved spacetime, there are two results known [2H20] concerning 
the non- negativity problem. If a solution of the EB system is given, in particular a four-dimensional 
manifold and a metric as a part of a solution, a result on non- negativity of solutions of the Boltzmann 
equation on suitable curved spacetimes implies a non-negativity result for the EB system. Thus in 
order to settle this question for the solutions obtained by Bancel and Choquet-Bruhat in [3] it suffices 
to examine whether the conditions assumed in [U [2D] hold in the framework of [3] . 

In Section 3.1, a result of Tadmon [20] is considered. Although [20] cites [3] it does not discuss 
the applicability of the results obtained to the solutions of that paper. The conditions in [5D] are 
quite complicated, but it is shown in Section 3.1 that they are satisfied in the framework of Bancel 
and Choquet-Bruhat. In Section 3.2 it is pointed out that the assumptions assuring non-negativity 
of solutions of the Boltzmann equation in the work of Bichteler j4i do not follow from those of [3] . 
It is shown, however, that non-negativity of the distribution function in the solutions of [3] can be 
proved using a method simpler than that of Tadmon. 

Section 4 contains some further general considerations on the EB system. The topics discussed 
include the advantages of a formulation using orthonormal frames, the issue of identifying physically 
relevant scattering kernels and how the non-negativity result can be made global. 

1.1 Notation 

In this part, we collect the notation which is used in this paper. Greek indices run from to 3, while 
Latin indices run from 1 to 3. The spatial variable x denotes a four- vector, while the momentum 
variable p denotes a three-dimensional vector, i.e., 

x={x°,x\x 2 ,x 3 ) 1 p^{p 1 ,p 2 ,p i ). 

The metric g a p has signature (—,+,+,+), g denotes its determinant, and the Minkowski metric is 
denoted by r\ a p. Throughout the paper, the speed of light c and mass of the particles are assumed to 
be unity and hence the momentum p a lies on a hypersurface defined by the equation g a pP a P^ = — lj 
which is called the mass shell. Due to the mass shell condition, p° and po are represented by p l as 
follows. 



p° = (gwp* + JigoiP*) 2 - goo{gijP i P j + 1) 

-ffoo \ v 



Po = -yigotP 1 ) 2 ~ 9oo(9ijP i P i + 1)- 
In some places, we use the Euclidean norm in M. d , i.e., 



(1.5) 



d 



W =2>T for per, 

and • will denote the usual inner product in M. d , i.e., 

d 

p ■ q = S~] ?>V f° r P,q ER d . 

We also collect the function spaces which are used in this paper and which were originally 
introduced in [3J. The spaces are basically Sobolev spaces weighted in the momentum variables. 
Two kinds of weight functions are used, (jp) N / 2 +\ k \ and e p , where iV" is a positive integer and \k\ is 
the number of derivatives which are taken with respect to momentum variables. For simplicity, we 
write hk(p) — (j>°) N ' 2+ > k \ and h(p) — e p , which can be understood as the case N — oo. 



1. Let wo be a domain of R 3 , i.e., x° = in R 4 , and wo = w o x R 3 - 

2. Let H^iujg) and H^ n(uiq) be the Sobolev spaces whose norms are defined as 

IMIh^o) = E H D * w llWo)> 

\k\<n 

\\f\\hM*o) = E n^(p)^, P /ii| 2(i0 ), 

\k\<n 

where k is a multi-index such that D% = D^D 1 ! in the latter case. 

3. Let n be a domain of R 4 and Cl = O x R 3 . 

4. Let Hfj,(il) and H^^n(Ci) be the Sobolev spaces whose norms are defined as 

IMIff^) = E W D xU\\h { n), 

\k\ <fi 
|fc|<M 

where fc is a multi-index such that D^ „ = D^D 1 ! in the latter case. 

5. g a (i{Q), dog a /3(0), and /(0) will denote the initial data for g a and / respectively. 

2 The Boltzmann equation in a curved spacetime 

In this section, we focus on the Boltzmann part of the EB system. We assume that a four-dimensional 
manifold and a metric are given, and consider the Boltzmann equation in this spacetime. Let 
x a £ I 4 , p" G R 4 , and p a satisfy g a pp a p^ = — 1. For a given metric gap, the Boltzmann equation 
reads 

Q(fJ)(x,p)= f f S(x,p 7 q,n)(f(p')f(q')-f(p)f(q))<m^dq. 
JR3 Js 2 V ' ~<10 

Here T l aj3 are the Christoffcl symbols defined by 

T M = ^9 aS {di9Sfi + dpgs-f - dggpj), 

and S is called the collision cross-section and is a non-negative function, 

S = Xga with a = a(g,d). (2.1) 

Q is a point of S 2 , thought of as the unit sphere in R 3 with standard measure dfl and polar angle 
8. The quantities a and 9 are called the scattering kernel and the scattering angle respectively. The 
scalar quantities A and g are defined by 



A = V-(p a + q a )(p a + q a ), g = V(p a -q a )(p a -q a )- (2.2) 



The variables p' a , q la , p a , and q a denote the momenta of two colliding particles, i.e., post-collisional 
momenta p' a and q' a , and pre-collisional momenta p a and q a , and they satisfy 

p' a + q' a =p a + q a , 

which expresses energy-momentum conservation. The collision operator Q can be written as Q = 
Q + — Q_ where 

Q+ = J J S(x,p, q, n)f(p')f(q r ) dn 1 -^- dq, 

Q- = f(p) If S(x,p,q,il)f(q)dQ^dq. 

Q + and Q_ are called the gain term and the loss term respectively. Note that Q+ and Q_ arc 
non-negative. 

In the case g a p = r) a p, the above Boltzmann equation reduces to the well-known special rela- 
tivistic Boltzmann equation. In this case, the Christoffel symbols vanish, \g\ — 1, and — qo = q°, so 
we obtain 

/ '" (>J - l f v 4 ,a(g,0)(f(p')f(q') - f(p)f(q)) dndq, 

Xg 



p a dx a 
where v^ is called the M0llcr velocity, 

V<f, 



p o q o- 



The general relativistic Einstein-Boltzmann system can be written in an alternative way which 
makes it look more similar to the special relativistic case and which is more convenient for some 
purposes. The idea is to introduce an orthonormal frame e" on spacetime with dual coframe d£. 
Here the indices from the beginning of the Greek alphabet like a are coordinate indices as before 
whereas the indices like [i from the middle of the Greek alphabet label the vectors of the frame. 
In formulating the EB system the metric is described by its components in a coordinate frame as 
before but the coordinates (x a ,p a ) used to parametrize the mass shell are replaced by coordinates 
(x a ,v u ) where v v = 0^p a . By abuse of notation we will write f(x a ,v u ) for the representation of 
the distribution function in the new coordinates. An advantage of the new coordinates is that the 
collision term no longer contains any explicit dependence on the metric and is identical in form to 
what it is in special relativity. On the other hand the transport part of the Boltzmann equation 
becomes more complicated and the Christoffel symbols are replaced by the components of the Levi- 
Civita connection in the chosen frame. Thus we obtain the equation 

"" e ^ - rFgapeZeUdseP + rf e e%)v\» ^ = Q(f, J). (2.3) 

If the EB system is to define a closed system for the quantities g a [j(x 1 ) and /(x 7 , v u ) then the 
orthonormal frame must be fixed in terms of the metric and the coordinates by some condition. This 
can be done in such a way that it depends only on the components g a p and does not contain any 
direct dependence on the coordinates. One way of doing this is to use the Gram-Schmidt process. 
This algorithm for producing an orthonormal frame from a general frame is best known in the case 
of positive definite metrics but it works just as well for pseudo-Riemannian metrics provided the 
starting frame is non-degenerate in the sense that the pull-back of the metric to the subspace spanned 
by any non-empty subset of the vectors of the frame is a non-degenerate quadratic form. This process 
can be applied to the coordinate frame d/dx a to get the frame e" Then the components e" are 
algebraic functions of the components g a p. In the case of the first two elements of the basis, for 



instance, we get 

e£ = #(-flbo)-* (2-4) 

e? = (-.9oo< + .goi^)((goo) 2 3ii - W.9oi) 2 )^ (2-5) 

It then follows that the components v a are also algebraic functions of the components g a p. The 
coefficients in the transport part of the Boltzmann equation are algebraic functions of the components 
g a p and their first order partial derivatives. 

2.1 The post-collisional momenta 

A parametrization of the post-collisional momenta p'^ and q'^ in general relativity has been obtained 
in [3] using an orthonormal frame {e^} associated to a point of the mass shell. Note that this is a 
different type of object from the orthonormal frame introduced above. In the frame considered in 
[3] p 1 * — e£v a , q* 1 — ejjuj" and v a + w a = (A, 0, 0, 0). The post-collisional momenta in this frame are 
represented by 

where fl £ S 2 . Consequently, p ,f * is given by 

p'^ = -(e^A + e^gcosO + gsm9(e2 cosy + e^ sinip)), 

where tt — (cos 9, sin 9 cos ip, sin 9 sin ip), and a similar expression is obtained for q 1 ^. It is a natural 
approach to choose an orthonormal frame and then use formulae from special relativity, but in some 
points it has an inconvenient aspect concerning the orthonormal frame. In this paper, we consider 
the Sobolev spaces of order greater than five. Hence, high order derivatives of some quantities are 
taken, and we need to estimate derivatives of the orthonormal frame chosen. If we follow the result 
of [3] , then derivatives of the post-collisional momenta are estimated as follows. 



\D k „ P p'\<(p°r^j:M^ 



P n 0\\r\ + \s\+l/2 



^ \ 2 \^\g 2 \^\\pxq\ 



g lM. 



where k is a multi- index such that D* = D^D^ the sum is over all the possible multi- indices 
ri, r 2 , s satisfying \r + s\ < \k\ and r\ + r 2 = r, and the term g^ k i> will be defined later. We 
can see that singularities appear in the denominator, which are certainly due to the choice of the 
orthonormal frame. This problem might be removed by taking another orthonormal frame at the 
singular points, but it would be quite complicated to find such an orthonormal frame that behaves 
well at those singularities. We instead use a different approach to parametrize the post-collisional 
momenta, where the above problem does not arise. 

We introduce a new parametrization of post-collisional momenta. Suppose that p a and q a are 
given, and consider the following four-vectors. 

n a =p a + q a and t a = (n i w i , -n Q w) for u £ S 2 , (2.6) 

Note that t a is orthogonal to n a . Then p la and q' a can be parametrized by 

P '" = p°- ^ (pg - gg) t°= P ° + 2Jg4t°, 

(2.7) 

q ~ q + t^ % ~ q z T^ t ' 



where we used t a n a — 0. It can be easily shown that they satisfy 

p' a + q' a =p a + q a and p' a p' a = q' a q' a = -1. 
We can use another parametrization for p' a and q' a , 



,,. _p a + q a g t a , a p a +q a g t a 



P ~ 2 + 2^F' q 2 2^' (2 ' 8) 



where g and t a are defined by (|2.2|) and (|2.6|) respectively. We can see that the parametrizations 
()2.7[) and (|2.8|) are very similar to those of the nonrelativistic case. In the nonrelativistic case, 
post-collisional momenta p' £ M 3 and q' £ R 3 are represented in two ways. 

p' = p — ((p — q) ■ u)u, q' = q + ((p — q) ■ w)w for oj £ S 2 

or 

, p + q \p-q\ , P + q \p - q\ , „ 2 

P = — ^— + — g— W, 9 = — g — W for W € 5 . 



Hence, (|2.7p and (|2.8|l can be thought of as natural generalizations from the nonrelativistic case. 

We remark that a well-known parametrization of p' a and q' a in the special relativistic case by 
Glassey and Strauss [9j is exactly the same as (I2.7[) . In the special relativistic case, g a p = 77^, we 
have 

t a t a - -(n ■ oj) 2 + (n ) 2 = (p° + q ) 2 - ((p + q) ■ uo) 2 , 

and by direct calculations, 

t a q a = -q°n ■ uj + n°U! ■ q = -q°(p ■ u) - q°(q ■ w) + p°(uJ ■ q) + q°(uj ■ q) 

= -p°q (p/p ■ u) + P°q°(q/q • w) = -p°q°((p - q) • w). 

Defining p = p/p° and g = q/q° we obtain the following expression. 

, 2p <f((p-q)-u)(p + q°) 

P P (pO + g 0)2 _ ((p + q) . w) 2 W ' 

This is the parametrization given in [9] . 

We close this subsection with a simple lemma, which shows that p a and —po are equivalent in a 
curved spacetime when its metric is close to the Minkowski metric. 

Lemma 2.1 Consider a momentum p a satisfying g a pp a p^ = — 1. Suppose that there exists a small 
e such that the metric g a /3 satisfies 

3 3 

\g a p-V a fi\<e and (1 - e)J2( xi f < dijX'X^ <(l + e)^pT) 2 

i=l i=l 

/or any i/iree dimensional vector X . Then we have 

|p|<(l+ei)min{p°,-jJo} and (1 - ei )p° < -p < (1 + £i)p° 

/or some small E\ > 0. 

Proof. By solving the equation g a pP a P^ — —1, p° and po can be expressed as functions of g^ and 
p 4 as in (ll.5[) . Hence, by direct calculations, 

bl 2 < (1 + e)<fapV < (1 + e)((.go,y) 2 ~ Soo (Si^V + 1)) = (1 + e)bo) 2 , 

-Po = -.9ooP° - goiP 1 < (1 + e)p° + e|p|, 
which imply that |p| < Cpo and |p| < Cp° for some constant C ~ 1. From the identity po = 
9oaP° + ffoiP 1 and the conditions on g a p, we get the desired equivalence between p° and — p$. D 



2.2 Derivatives of the post-collisional momenta 

In this paper, we use (|2.7[) for p' a and q' a . Since the parametrization (|2.7[) depends on x, p, and q, 
derivatives of the post-collisional momenta are not trivial to compute, but we can see that derivatives 
of p' a and q' a depend only on t a . To calculate derivatives of p' a and q' a , we use the notation <?( m ) 
of Bancel and Choquet-Bruhat from [3] . 

A function h is said to be a g( m ' function, if h is any linear combination of products, FT . D k J g^, 
with ^ \kj\ < m and coefficients in the algebra of bounded functions on fi x M 3 x R 3 generated by 
p a /p°, p a /p , q a /q°, and q a /q . 

Lemma 2.2 Let t a be a four-vector defined by 
multi-index k = k + k, we have 



for some w G S . Then, for D k — D k D k with 



D k 






2+2|fc 



= (P°r lkl E 



(=0 



(p0)I(g0) 3 + 2l fcM (ISI) 

{t a t a Y+\ k \ 9l ' 



where g\ are some g" k \> functions. 

Proof. We first note that p-derivatives of p a and p a are g^ functions as follows, 

dp pk p a dp 1 

Qpk 



Pk P 

' = 9ak, 

Po Po 



dp k 



>k- 



dp a dp 



dp 1 



—L- — TTTSlOa + -TT-rgia — <7/3fcffOa + 9ka , 

op K ap K op K 



p" 
Po 



while their ^^-derivatives are p°g ( - 1 ' functions, 

dp pPp 1 dpp dp 

dx^ = ~^o~ da9 ^ d^ = dx^ 9 ^ +P" d "^ 

Hence, we can deduce that 



'^Po~ 



d a g~/6goi3 +P 1 d a g t 3 1 . 



dgb 



9 



(m) 



and 



dg (r 



Qpk p0 Q x a 



_ (m+l)_ 



To prove the lemma, we use induction on |fc|. Since n Q = p a + q a , the four- vector t a is of the form 

t a = {n.uj 1 , -n uj) = nP{gipw % , -g puj) = p°gQ 0) + q°g[ 0) 
for some ^g and g\ , which are g^ functions. Hence, we obtain for \k\ = 0, 



W*, ^(P°)V) 



tptP 



' Q = E 



0W„0\3— i 



t t a 



-g 



(0) 



n 0W„(n2-i JO) 



We use t a t a = T,Uo(pT(l ) 9i to show that for l fc l + L 

2+2|fc| 



D„cD k 



2^ 

tpt? 



= D X « 



f _p r i*i v ( P °r( g °) 3+2|fc| - ( | £ |) 



/=o 



(t a t a y+\ k \ 



<p°r* 2+ £' m^ e^ 



8=0 

2+2|fc| 



(*/^) 



1+lfcl 



/V,0\if„0\3+2|fc|-i 

(p°)- |fc| E M ( LL, a-[/-/ 



,(l*l) 



i=0 



(i^) 2 +l fc l 



Since D xa [tptP] = Ei=o(P )*(9°) 2_i ffi > we obtain 



Similarly, 



D^D K 



2^V 



. , 4 ±l fe i ^0)i^0-)5+2|fc|-i 



^-ifci y: 



j=0 



(^^) 2 +l fc l 



-.9 



(I*1+i) 



DpiD* 



2^ 



= A, 



2+2|fc| , si( 0\3+2\k\-i T 



= (p )- 1 *"- 1 E 



(t Q i«)i+ifei 

2 + 2 l fc l ( 0W„0\3H-2|fc|— t 



(p U )%°) 



foi^+W 



,(|fc|) 



2+2|fc| 



(p°)- |fc| g (P (^) 2+ i fc i ^[^]ft 



(m) 



i=0 



Since # pJ [t^] = (p )- 1 ^(p )'^) 2 -^ ', we obtain 



D xa D k 



2^' 

^f 3 



(P' 



o^— |fe|— x \p iP)w) (|fc|) 

(t^)2+|fc| « ' 



i=0 



and this completes the proof of the lemma. 



□ 



To estimate derivatives of the post-collisional momenta, we need to obtain a lower bound for the 
quantity t a t a , which appears in the denominator in the expressions coming from Lemma 2.2. Note 
that t a t a is non-negative, because t a n a — and n a = p a + q a is a timelike vector. 



Lemma 2.3 Let t a be a four-vector defined by (|2.6j) for some uj 6 S 2 . Suppose that there exists a 
small e such that the metric g a p satisfies 



\9ap-V*p\<£ and (1 - e) E(*T < 9ij^X j < (1 + e) ^(X*) 2 



for any three dimensional vector X , then we have the following lower bound. 

M 2 + (go) 2 ' 
3p qo 



t a t a > 2( 50l o/r - 9oo(9ijOJ 1 ^) 2 + 



Proof. The proof is a direct calculation. Since t a = (riju/, —uquj), we have 



i\2 



t a t a = (n ) (gij^ui 3 ) - 2(n i uj l )no(g oi uj l ) + g 00 (niUJ l ) 

= (n Q f{g ii u i uJ i ) +9oo((fi'ynV) 2 + 2( 3lJ nV)n ( 5oi w i ) + (n^fooiW*) 2 ) 

- 2{(g ij n % u) 3 )n (g 0i u] 1 ) + n°n Q (g 0i uj') 2 ), 

where we used riiUi 1 = g^n 1 ^ + n°(goiLu l ). There are six terms on the right hand side of 
first two of these terms can be estimated as follows using the Cauchy-Schwarz inequality. 



(2.9) 
. The 



(no) (gijUJ 1 ^) + goo(gijn l uj ] ) > (no) (gij^u 1 ) + g 00 (g lJ n t n : ')(g ij u} l u} : '), 



We continue the estimation as follows 

(no) 2 {gijU % u> j ) + goo^jrtuj 3 ) 2 

> (g^uj'uj^dpo) 2 + (g ) 2 + 2p <2o + ffoo(ff«pV) + ffoo(5«gV) + 2g m {g l0 p L q 3 )) 

= (flfywVjCCgoiP*) 2 + (5«9 4 ) 2 - 2goo + 2p 9o + 2g 00 (g ij p i q j )), (2.10) 



using 



po = -\j{goiP 1 ) 2 - 9oo((j9ijP i P j ) + l) 



and the analogous formula for go- The fourth and the sixth terms on the right hand side of (|2.9|) 
are calculated as follows. 

(ffo 4 w l ) 2 (5oo(n ) 2 -2n n ) 

= (,9o^) 2 (5oo(;p ) 2 + ffoo(g°) 2 + 2gooP°q° - 2p°p - 2p°q - 2q°p - 2g%) 

= (goi^fi-gooip ) 2 - goo(q ) 2 - 2WV - 2 P °(g 0i n i ) - 2q°(g 0i n% 

where we usedpo = gooP° + (goiP z ) and the analogous formula for qg. To proceed with the calculation, 
we use the formula 

p° = ( goiP 1 + \ {goiP i ) 2 ~ goo(gt 3 p l p J + 1) 

-.9oo V v 

and the corresponding formula for q°. After a long calculation we obtain 

( 50l ^) 2 ( 5 oo(« ) 2 -2n n ) 

= ( 90l wf(( 5l3 pV) + (.9 4J gy ) + 2) 

+ 2(g 0i u i f J{gmP % ) 2 ~ goo((g^P l P f ) + l)J{g 0l q 1 ) 2 - .9oo((s*j9V) + 1) 

—5oo v v 

- 2{g 0i w i f-^—{g 0i p i ){g fH q i ). (2.11) 

-300 

The third and the fifth terms on the right hand side of (|2.9[) are estimated as follows. 
(goiV l )(gijn l u j )(2g on° - 2n ) = -2(goiUJ l )(gijn l u] j )(g in l ) 

> -2\gaiW % I J gijrini y g ij oj i uj^\goin l 

> -(W) 2 (s«nV) - (ftiwVjCswn*) 2 , (2.12) 

where we used the Cauchy-Schwarz inequality. Hence all terms on the right hand side of (|2.9p have 
been estimated by (|2.10p - (|2.12l) . Moreover, since n a — p a + q a , the first two terms on the right 
hand side of (|2.10[) are partially cancelled by the second term on the right hand side of (12.12[) . and 
similarly the first two terms on the right hand side of (|2.11j) are partially cancelled by the first term 
on the right hand side of (|2.12|) as follows. 



and 



fewV)(( ff oiP i ) 2 + (g 0t q 1 ) 2 ~ (ffo^) 2 ) = -2( 5 ^V)(.9o^)W), (2-13) 

(soiwftfftjP'V) + G?ij«V) - (ffynV)) = -2(ffoi« i ) 2 (g ii i»V) 

> -2(« ?0l ^ l ) 2 V /^Vy / 5^9V. (2.14) 
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As a consequence, (|2"1?)) is estimated by (|2.13|l - (|2.14|) and the other terms of (|2.10p - (|2.1ip . i.e., 
t a t a > (g l3 uj l uj : >)(-'2gQ + 2p q + 2g 00 (g ij p i q j ) - 2{g oi p l ){g Qi q 1 )) 
+ (ffo*^) 2 (-2yfgijP i P j yfgij<l i <l j + 2 J 
+ 2{g 0i u i ) 2 J{g 0i p i ) 2 - goo((5yPV) + l)J{g Qi q 1 ) 2 - Soo((3uW) + 1) 

— Onn v v 



-300 

2( 50l u/) 2 (go i p t )(go 1 q l ), 

-goo 



and we rewrite it as follows. 



M Q > -2(flfywV>oo + 2(.go,:w 1 ) 2 

+ 2(g lJ u l uj : >){poqo + g Q0 (g l:j p l q : ') - (goiP l )(gmq 1 )) 

+ 2(g 0i uj i ) 2 JigaiP*) 2 ~ 3oo(toV) + 1) J {go l q 1 ) 2 - Soo(taV) + 1) 

— Onn v V 



-.900 



- 2( 50 iW i ) 2 — — (-gooyJgijP'piyJgijq'q' + (W) (#>*«*) J • (2.15) 

Note that the first line of (|2. 15|) is strictly positive, and the sum of the third and the fourth lines is 
non-negative because of the following calculation. 



-gooy gtjP^ y gijq i q j + (goiP l )(goiq' t 

= {gaof {gi]P l p 1 ){gijq l q J ) + {goip 1 ) 2 '{gm 1 ) 2 - 2goo(j9oiP t )(9oi9 t )ygijP i p j y9ijq i q j 
< (fldo) 2 (ftjpV)(fti«V) + G?0ii» i ) a (9di« < ) a -goo{{gotP l ) 2 (gijq l q j ) + (W) 2 teV))- 

We now have only the second line of (|2.15[) . which is explicitly calculated as follows. Note that the 
coefficient gijofuj 3 is strictly positive, and we have 

poqa + gooigijpW) - (gotpKgmq 1 ) 

>Poqo- [ -goQ\JgijP i p j \Jgijq i q j + \goiP l \\goiq l 



{poqo) 2 - \-gm\J gi]P l p 3 \J gi^qJ + \goiP t \\goiq l 

iP'Wgoiq 

1 ( / / — / — 

> 



pogo + i^-goo\/g7jP T p 1 ^/gij¥¥ + \goi. 
3 ( (Poqo) 2 - f — ffoo-y/ j?ypV y 9ijq i q j + \goiP l \\goiq i \ J J , (2.16) 



where we used Lemma 2.1. Note that 

(poqo) 2 = (go i p l f{go l q l f 

- goo{goiP l ) 2 {gijq l q j ) - goo(goip') 2 - goo{goiq % ) 2 {gaP % p j ) - goo(goiq l 
+ {goo) 2 {{gijP l p i ){g l] q l q j ) + gijPY + g ij q i q j + 1), 
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and 



-9ooygijP i p j y9ijq i q j + \gmp l \\gmq l 
= (W 2 (5yPV Xft^y ) - 2 9oa\goiP t \\9oiq' l \y gtjP^ 3 y gijq l q 3 + (goiP t ) 2 {go l q 1 ) 2 



< (9oo) 2 (gijP i p i )(gijq i 9') -goo{goiP i f{gij<i i <i') - 9oo(goiq i ) 2 (9ijP i p i ) + (goiP l ) 2 (gotq 1 ) 2 - 

Then, (I2.16[) is estimated as follows. 

poqo + goo{g lJ p l q r ) - (gotP l )(gotq l ) 



> 



3poqo 



(-gaoigoiP 1 ) -goo(gozq l ) + (goo) {g^p 3 +gzjq t q J + l)) 



>^((po) 2 + (qo) 2 ), 



(2.17) 



where we used the explicit formula (jl.51) for po and qo- We apply the estimate (|2.1T[) to (|2.15p . and 
obtain the desired lower bound. D 



Remark 2.1 In the special relativistic case, we have a more refined result, 



t a t a > 2 + \p x uj\ 2 + \q x lo\ 2 + 



1 + H 2 



VT+Wy/T+W 

which was obtained in J§j/. This inequality implies that t a t a > p° /q in the special relativistic case, 
and Lemma 2.3 shows that a similar inequality holds in a curved spacetime, i.e., 

t a t a >C Po /qo or t a t a >Cp°/q° 



(2.18) 



for some positive constant C . 



We apply the inequality (|2. 18|) to the result of Lemma 2.2, and obtain the following estimate. 
Under the same assumptions as in Lemma 2.2, we have 



D k 



2^e 
. tpt? 



<(p u ) 



2+2|fc| 
"■-I :i'l. V^ C„0\J+|fc|/"„0\4+3|fc|-J 



E (p y +lkl (q°) 



>=o 



,(|fe|) 



<C(p°) 1+ W(q°) 4+ W 



,(ifci: 



for some g" k V and a positive constant C. Since p' a and q' a are parametrized by (J2.7I) and it can 
be easily shown that D k p a = (p°) l ~' k ' g" k '> , we obtain the following result on derivatives of the 
post-collisional momenta. 

Lemma 2.4 Let t a be a four-vector defined by (|2.6p for some u G S 2 , and p' a and q' a be two 
post-collisional momenta parametrized by (J2.7I) . Suppose that there exists a small e such that the 
metric g a p satisfies 

3 3 

\g a p-il a p\<e and (1 - e) £pT) 2 < g^X' < (1 + e) ^(X 1 ) 2 

i=l z=l 

for any three dimensional vector X . Then, derivatives of p' a and q' a can be estimated as follows. 

,(1*1) 



\D k p' a \ + \D k q' a \ < C(p°) 1+ W(qy +3 W 
for some g(\ k \> and D k = D k D k with multi-index k = k + k. 



,9 V 
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2.3 The /i — N regularity of the collision operator 

In [3J, the authors introduced a regularity property of Q called \i — N regularity, and showed that 
the Cauchy problem for the EB system is well-posed when Q satisfies it. The collision operator Q 
is said to satisfy [i — N regularity if there exists a constant C satisfying for each ui t , 



flQ(/,/)(t) 



<c\\f(t)\\%. 



Hp,N{&t) 



n,iv(&t) ' 



To obtain jx — N regularity, the collision cross-section S should satisfy suitable conditions, and by 
following the calculations of [3], we can see that the following inequality also holds under the same 
conditions on S. 

' ' ' <C||/(i)||^(* t )||<?(*)lk,^), (2-19) 



-Q(Lg)(t) 



where Q{f,g) is defined by 

Q(f,g) = 



S(x, P ,q,n)(f(p')g(q')-f( P )g(q) 

\ 



dQ dq. 

-9o 



In the previous section, we obtained a new estimate of D k p' a and D k q' a in Lemma 2.4. This 
leads to a corresponding new set of conditions on S which imply the \x — N regularity of Q. Note 
that high order derivatives D k p [f(x,p')] are linear combinations of the following quantities. 

(D\f)(x,p')(D^p')---(D^p') with i 1 + .--+i«=|fc|, (2.20) 

where D l denotes some differential operator D r x satisfying \r\ = i. Lemma 2.4 gives the estimate 



\Dl p [f(x,p')]\ < C^|P7)(^p')l(p ) i+|I| (« ) 4i+3|S 



,(|fc|) 



where the sum is over all the possible i satisfying (|2.20[) . Consider the following quantity which 
arises from estimation of the gain term Q + . Let the multi-indices k, r, and s satisfy k + r + s = I 
with \l\ < fi. 



1 



1.91 



t) n I J {Dl p S)D^ p [f(p')}Dl p [f{ q ')}dn ] ^-dq 



<cJ2[ \Dl P s\\D*f( P ')\(p y+^(q°r^ 



7 (l?l) 



x\&f( q ')\(p° 



P\J+\S\f q 0\43+3\S\ 



<cY, K 2 (p)hW)\D l !(p't 



h](q')\nif{q>)\ 



,(|5|) 



,(|5|) 



dn- 



Po9o 



■ dq 



'dn^dq 

Poqa 



x // |^, J>1 s'| 2 ^(p)/ it - 2 (p')^7 2 (e')(p ) 2(i+i)+2|r+s| (g ) 8(i+i)H - 6k+s| dn — 



i 



Pbq Poqo 



dq, 



and if the second integral is bounded, then we multiply the first integral by the weight function 
hf(p) and integrate it over x and p with ^^ dpdq = dp' dq' , to obtain the inequality (|2.19[) when 
H > 5. Hence, the boundedness of the second integral should be the condition on S for the /i — N 
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regularity of Q. The second integral is estimated as follows. In the case where N is finite, 

J Po%PoQo 

<C \D k S\ 2 ^X r (p ) 2 lH-l+2|r+.|-2/^) ) 14|r+.|-l d n d g 



<C | J D^5| 2 (p°) 2 l'l(p°) 2 l r+s l+ 2 l'" +s l- 2 ( ? ) 14 l , "+ s l- 1 ( iO^ 

< C ff\D k ^ p S\ 2 (p ) 2 ^(P°f r+S ^ 2 (q ) lMr+shl dndq, (2.21) 
where we used — p < Cp' Q qQ and p° < Cp'°q'°. For the case N = oo, 

\D k x v S\ 2 h 2 (p)h- 2 (p')h- 2 (q')(p°) 2 ^^~\(qy^^~Un ^ dq 

Po%poqo 

< O If \D k xp S\ 2 h- 2 (q)(p°) i \ r+s \- 2 {q G ) li \ r+s \- 1 dfl dq, (2.22) 



where we used p'° + q'° = p° + q a and — p < Cp' a q' Q . The boundedness of the integrals (j2.21|) 
and (J2.22I) are the conditions on S for the gain term. The arguments for the loss term are much 
simpler, and it is easily shown that the condition for the gain term implies that for the loss term. 
Consequently, we obtain the following lemma. 

Lemma 2.5 [3] Let /j, > 5, and suppose that the collision cross-section S satisfies the following 
conditions for multi-indices k + r = I satisfying \l\ < \i. 



oc. 



l^^| 2 (g°) 14|rhl ^^<C(p°) 2 - 4|r| - 2|fe| for N< 
\Dl p S\ 2 e- 2 <'\q ) u ^- 1 (mdq<C(j> ) 2 - i ^ for N = 

Then, the collision operator Q satisfies /j, — N regularity. 

Remark 2.2 Note that the condition on S for finite N is stronger than that for the case N = oo, 
and that the former condition implies the latter one. 

JJ\D k x \ p S\ 2 e- 2 ^{qy^-'d^ldq 

< f[\Dl p S\ 2 (q°) u ^- 1 dndq < C(p°) 2 - 4 M-2|fc| < C(p ) 2 - 4 ^. 



In other words, if the collision cross-section S satisfies the condition for finite N, then it also satisfies 
the condition for the case N = oo. 

3 The non-negativity problem for the Einstein-Boltzmann 
system 

We now consider the EB system and the question of the non-negativity of its solutions for non- 
negative initial data. The Cauchy problem for the EB system was first studied by Bancel and 
Choquct-Bruhat [21 IS]- We apply Lemma 2.5 to the results of [UG3 to obtain the following existence 
theorem. 
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Theorem 3.1 [H|3] Suppose that Cauchy data are prescribed as 

j a/) (0)e^ +1 (w o ), <9 O 5a/3(0) e Hn(u ), /(0) € JT m ,jv(wo), 

swc/i t/iat 

|g Q/3 (0) - ?7q/3 | < e - <5, <5 > 0. 

If 1^ > 5, JV > 6, and i/ie collision cross-section S satisfies the conditions of Lemma 2.5, then there 
exist a domain £1 in M. 4 , which admits luq as a Cauchy surface, and a function f on Cl such that 

1- g a /3 e H/i+iipi) and f e H^^Cl). 

2. g a/3 satisfies \g a/3 - r) a p\ <e on O. 

5. <?a^ cmd / satisfy the EB system. 

4- 9ap an d f induce the prescribed Cauchy data on uj and uj respectively. 

This solution is unique in Q and depends continuously on the Cauchy data. 

In this section we consider the non- negativity of the distribution function constructed in the above 
theorem. Note that we already have a solution (/, g a p, fi) of the EB system, while the non-negativity 
problem concerns only the distribution function /. Hence, it is enough to consider the Boltzmann 
equation on a given curved spacetime and show non-negativity in this case. Conditions for the non- 
negativity of solutions of the Boltzmann equation on a given spacetime have been given by Bichteler 
[4] and Tadmon |20j . We will investigate to what extent the arguments for non-negativity given in 
[4j [20] apply to the solutions of the EB system constructed in [3] . 

Theorem 3.1 holds for any weighted spaces with N > 6, but we will only consider the case 
N = oo, i.e., the exponentially bounded case. Firstly, any functions that have polynomial decay 
rates can be approximated by functions having exponential decay rates. Hence, non-negativity of 
the case N = oo implies the corresponding result in the case of finite N by suitable approximations. 
Secondly, as we checked in Remark 2.2, a collision cross-section which satisfies the condition for the 
case of finite N also satisfies the condition for the case N = oo. So, any property that holds under 
the condition for the case N — oo holds under the condition for the case of finite N. 

3.1 The result of Tadmon 

The original method for the non-negativity problem introduced by Lu and Zhang 10J has recently 
been applied to the general relativistic case by Tadmon 20 . A curved spacetime was assumed to be 
given, and mild solutions of the Boltzmann equation were considered, which are defined as follows. 
Consider the vector field (p a , -r« pV) as in [4], and parametrize its integral curves by x° — s on 
the mass hyperboloid. Let (X(s),P(s)) £ M 3 x R 3 denote an integral curve with X z (0) — x l and 
P(0) — p parametrized by s. This has the physical interpretation of a particle path. Along this 
curve, the Boltzmann equation is written as follows. 

f{t, X(t), P(t)) = /(0) + f K(f)(s, X(s),P(s)) ds, (3.1) 

Jo 

where /(0) is an initial datum evaluated at (x % ,p), and K denotes 

K(f)(x,p) = ^Q(f,f)(x,p). 

Mild solutions are now defined as follows: a function f is called a mild solution of the Boltzmann 
equation with measurable initial value /(0), if f is measurable, K(f) is Lj- oc (R + ) along the integral 
curves, and p. II) holds. 
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For the rest of this section, x will denote three-dimensional vector consisting of the spatial 
components x % where 

x a = (cc ,^) = (t,x) £l+x M 3 , 

and the distribution function will be written as f(t,x,p) instead of f(x,p). The main theorem of 
[20] can be stated as follows. 

Theorem 3.2 Let f be a mild solution of the Boltzmann equation with a non-negative initial datum 
/(0). Assume the following conditions. 

S(t,x,p,q,il) _ S{t,x,q,p,n) 

p'qo q u Po 

S(t,x,p',q',Q.) _ d(p,q) S(t,x,p,q,Cl) 

W P'Wo ~ d(p>,q>) p Q q 

r r q o 

(hi) ess sup // 9 \f(t,x,q)\\g\? dfldq =: a(t) £ Ll c (R+), 

x,p JJm^xS 2 (P ) 10 

(iv) esssup^M4Z(^). :/3(i)GLL(M+) , 

where I is defined by 

l(s,X(s),P(s)) := d^isp-^Xis^Pis)) 

with J, the Jacobian for 

d(x,p) 



dx dp - 



d(X(s),P(s)) 



dX(s) dP(s) =: J(X(s), P(s)) dX(s) dP(s). 



Then, f is non-negative. 



Remark 3.1 The conditions on a(t) and j3(t) should be replaced by a, (3 £ L* oc (I) for some finite 
interval I, because the solutions in Theorem 3.1 are local in time. Hence, it is enough to show the 
boundedness of a and ft on I . 

In this part, we will show that, with an extra assumption, the solution of Theorem 3.1 satisfies 
the conditions (i)-(iv). Since, as explained below, the extra condition can be arranged by making a 
coordinate change, non-negativity for the EB case is ensured. We remark that the strong solutions 
in Theorem 3.1 satisfy the definition of mild solutions. 

Non- negativity of /. We consider the conditions (i)-(iv) separately in the following. 

(i)— (ii) The collision cross-section S is given by (|2.ip . which implies 

S(t,x,p,q,n) = S(t,x,q,p,n) = S(t,X,p',q',il), 

and it is well known that the Jacobian in (ii) is given by 

dip, g) _ pm 
9ip',q') p' q' ' 

Hence, the conditions (i)-(ii) can be written as po/p° — qo/q° = Pq/p'° for any p and q. However, 
this leads to 

g 0l = 0, (3.2) 



1G 



which is not assumed in Theorem 3.1 and is the extra condition referred to above. 



(iii) For the metric g a p given in Theorem 3.1, its determinant \g\ is bounded, and go and q° are 
equivalent. Therefore, (iii) is estimated as follows. 



e 2q °\f(t,x,q)\ 2 dttdq 



a(t)<Cj-^JJ S\f(t,x,q)\d£ldq 



<C 



(P ) 2 



S 2 e- 2q ° dfldq 



where e q is the weight function for N = oo case. Since / 6 H^ t N(Q) with /i > 5, the second integral 
above is bounded. The condition on S given in Lemma 2.5 implies the boundedness of the first 
integral, and this shows that the condition (iii) holds in Theorem 3.1. 

(iv) As for the last condition, we recall that the integral curve (X(s),P(s)) is defined as follows. 
Let X(s) and P(s) denote for each x and p, 



and they satisfy 



X(s) = X(s, 0, x,p), P(s) = P(s, 0, x,p), 



P a (<s)pP(s) 

d s p\s) = -r aP po (a) . pl (°) = p 1 



(3.3) 



where the Christoffel symbols are evaluated at X a (s) — (s,X(s)). We first consider the quantity 

p°(s)] = (d x *p°)(s)(d s x<*( s )) + (d p .p a )( s )(d s p*(s)) 



pP(s)Pi(s) 
2P (s) 



{d a gi3~ 



P a {s) Pi(s) j P a (s)P^(s) 
P°{s) Po(s) ap P°(s) ' 



where we used d x c,p° = —p^p 1 d a gp 1 /{2p ) and d p kp° = —pk/po as in the proof of Lemma 2.2 with 
the equations (|3.3[) . Since g a p <E iJ p +i(0) with \x > 5, d a g^ and T l a g are bounded. Lemma 2.1 
then gives the estimate 



d* 



P°(s) <CP°(s). 



(3.4) 



Consider now d(X(s), P(s))/d(x,p). We need the following calculations. Differentiate the first 
equation of (|3.3[) with respect to x 3 to obtain 



QQ X i (s) 9^(3) P i ( S ) 

_ d x] P l (s) P\s) 



P°(s) 
P^(s)Pf(s) 

P0( S ) (P0( 8 ))2 V 2P (s) 

and as a consequence we obtain the estimate 



(d a g 01 )(d xl X a ( S )) 



Pj(s) 
Po(s) 



(d^P'is)) , 



Similarly we have 



d s \d x X(s)\ < C\d x X(s)\ + C—-\d x P(s)\. 



d s \d p x( s )\ < c\d p x(s)\ + c„\d p P(s)\. 



(3.5) 
(3.6) 
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Differentiate the second equation of (I3.3[) with respect to x 3 to obtain 



Q ^ (P°(s)) 2 a - J 



P°(a) 
P°( S ) 



a/3 



P°(s) 



P a ( S ) 



By similar arguments, we have 

d s \d x P(s)\ < CP°(s)\d x X{s)\ + C\d x P(s)\, 

and similarly again 

d s \d p P(s)\ < CP°(s)\d p X(s)\ + C\d p P(s)\. 

Due to (|3.4j) . the inequalities (J33J-((3?H]) can be combined as follows. 



(3.7) 
(3.8) 



|P°(*)0 x jr(s)| + |0 B P(s)|] < c(|P°(a)5 B Jr(a)| + |0 B P(*)| 
\P°(s)d p X(s)\ + \d p P(s)\] < c(\P°(s)d p X(s)\ + \d p P(a)\ 



We now use Gronwall's lemma with the conditions X(0) = x and P(0) — p to obtain 
\P°(s)d x X(s)\ + \d x P(s)\ < c(\P°(0)d x X(0)\ + \d x P(0)\) < Cp°, 
\P°(s)d p X(s)\ + \d p P(s)\ < C(|P°(0)5 P X(0)| + \d p P(0)\) < C, 

which imply the estimates 

\d x X( S )\< Cl ^- y \d p X( S )\< Cl ±- y 
\d x P(s)\<Cp°, \d p P(s)\<C. 



(3.9) 



If we apply (|3.9[) to (|3.5[) - (|3.8[) again, we can see that s-derivatives of the above quantities satisfy 
the same estimates. 



d s \d x X(s)\ < C-^- y d s \d p X(s)\ < Cj^, 



(3.10) 



d s \d x P{s)\<Cp°, d s \d p P(s)\<C 
On the other hand, we can consider the following integral curves. Let Y(s) and V(s) denote for 



each y and v, 
and they satisfy 



Y{s) = Y{s, t, y, v), V(s) = V{s, t, y, v), 



VHs) 



v°{ s y 

where the Christoffel symbols are evaluated at Y a (s) — (s,Y(s)). Then, by the same arguments as 
above, we obtain 



\V°(s)d y Y( S )\ + \dyV(a)\] < c(\V°(s)d v Y(s)\ + \d y V(s)\ 
\V°(s)d v Y(s)\ + \d v V(s)\\ < c(\V°(s)d v Y(s)\ + \d v V(s)\ 
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and use Gronwall's lemma with Y(t) = y and V(t) = v, 

\V°( S )d y Y( s )\ + \d v V(s)\ < c(\V°(t)d y Y(t)\ + \d v V(t)\) < Cv°, 
\V°(s)d v Y(s)\ + \d v V(s)\ < c(\V°(t)d v Y(t)\ + \d v V(t)\) < C, 

and as a consequence the following estimates are obtained. 

1 



\d y Y(s)\ < C- 



\d v Y(s)\ <C- 



V°(s)' ,v yn ~ V°(s)' 
\d y V(s)\ < Cv\ \d v V{s)\<C. 

We now consider the quantity 

/( s ,x( s ),p( s )) = a s [p ( s )]j- 1 (x( s ),p( s )) + p°( s )a s [j- 1 (x( s ),F( s )) 

from which we have for y — X{t) and v = P(t), 

l{t, y, v)J(y, v) = d s fp°(s)l + v°J(y, v)d s \j^(X(s), P(s)) 

J s—t L 

By applying (13.41) to the first quantity above, we obtain 



d s 

Note that J(y, v) is written as 

J(y,v) = J(X(a),P(a)) a=t 



P°(s) 



<Cv°. 



d(x,p) 



d(X(s),P(s)) 



d(Y(s),V(s)) 



d(y,v) 



where we used 

x = Y(0,t,X(t),P(t)) = Y(0,t,y,v)=Y(0) and similarly p = V(0). 

Thanks to the multilinearity of determinant, J(y,v) is estimated by (j3.11[) as follows. 

„o\3 



J(y,v)<C 



v° \ 3 =r( — 

v°( s )J s=0 U° 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



On the other hand, as for J _1 (X(s), P(s)), we have the same estimates for d x X(s) and d s d x X(s), 
d p X(s) and d s d p X(s), and so on, as in (I3.9[) and (|3.10l) . Thanks to the multilinearity again, the 
following estimate is obtained. 



ds 



J- l (X(s),P(s)) 



<C 



P°(s) 



C 



(3.15) 



We combine ([3~T3l> - (|3~T5l) and (f3TT2|) to obtain 

\l(t iy ,v)\J( yi v)<Cv°. 

This gives boundedness of /3 on /, and therefore it is proved that the condition (iv) holds in Theorem 
3.1. 
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It remains to discuss the extra condition g$i = which is required to ensure that the conditions of 
Tadmon's theorem are satisfied. We are concerned here with a solution on a local region, where the 
components of the metric in a certain coordinate system are close to those of the Minkowski metric. 
Let the original coordinates be (t, x 1 ). Now choose new coordinates (t,x l ) with the properties that 
x l agrees with x % for t = and the x l are constant along the integral curves of the unit normal 
vector field to the family of hypersurfaces of constant t. Then the components of the metric in the 
new coordinate system satisfy the desired condition. The new coordinates might not be defined 
on exactly the same region as the old ones so that non-negativity is obtained on a slightly smaller 
domain. Given the fact that the non-negativity theorem being discussed here is local in nature, this 
is not a major disadvantage. 

3.2 The result of Bichteler and a new approach 

In Bichtcler's paper [J], it was proved that under certain hypotheses a local solution of the Boltzmann 
equation exists on a given curved spacetime, and that it is non-negative. A four-dimensional manifold 
M and a Lorentzian metric g a p were assumed to be given, and the Boltzmann equation was written 
in the form 

Cxf = IIJ W(12 -> 34) 5(1 + 2 - 3 - 4) (/(3)/(4) - /(l)/(2)) 32 d'i d4, (3.16) 

where 1, 2, 3, and 4 stand for p a , q a , p' a , and q' a respectively. The left hand side of (|3.16[) means 
that the derivative of / along the integral curves of the vector field X = (p a , — T^ p^p 7 ) exists at 
almost all points of the domain of /. The transition rate W is given as 

W(12-»34) = fccr(12^34), 

where k is some kinematical factor, and a can be regarded as the same quantity as the scattering 
kernel defined in (|2.ip . The volume forms d2, d3, and dA are same as in our case. 

d2 = ^dq, 

-<7o 

and d3 and dA are similarly defined. The main theorem of [4] is the following. 

Theorem 3.3 [4] Let /(0) be a measurable junction on Cjq, and suppose that there exists a contin- 
uous timelike vector field fi a (x) on ujq such that 

f(x,p)<Ce^ x '> pa on wo, 

and the scattering kernel satisfies the following property. 

cr(12 -> 34) 5(1 + 2 - 3 - 4) di dA < const. (3.17) 



Then, there exists a solution f to the Boltzmann equation (|3.16j) on a domain Cl, which is again 
exponentially bounded. Moreover, if /(0) > ; then also f > 0. 

The condition (|3.17[) can be stated in our notation as 

a(g, 0) dQ, < const. (3.18) 



The corresponding conditions on the scattering kernel for the EB case are given in Lemma 2.5 
with the relation (I2.1J) between the scattering kernel and the collision cross-section. Although the 
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conditions in Lemma 2.5 are much more complicated than (I3.18|) . we can see that they do not imply 
(|3.18p . Hence Theorem 3.2 cannot be applied directly to prove Theorem 3.1. 

If there exists a solution of the Boltzmann equation with non-negative initial data which becomes 
negative at some point in phase space then it is natural to follow the particle path through that 
point backwards in time. This curve must meet the Cauchy surface and thus there is at least one 
point on it where / is non- negative. Hence there exists a point (£»,x*,p*) on the curve where it is 
zero and immediately afterwards negative. At this point the derivative of / in the direction of the 
particle path is non-negative because the loss term vanishes. It might however be zero. This type 
of consideration plays a role in the proof of positivity in [3] 

A small modification of this idea can be used to give a relatively simple proof of non-negativity 
for the EB system. We first introduce the following notation, 

d P^ d , p a p li d 



dt p° dx l a0 p° dp 1 ' 

and write the Boltzmann equation as Cf = {p a ) 1 Q(f, /)■ We then modify the Boltzmann equation 
by adding a small quantity on the right hand side. 

^fn = ^Q(JvJv) + Ve-^\ (3.19) 

where r\ > is a small parameter. The above equation is understood to be defined on Cl with a 
metric g a p constructed as in Theorem 3.1. Since the quantity rje~' p ' is sufficiently smooth and 
square integrable on O, existence of solutions is easily proved by the same argument of [5] using the 
following lemma. 

Lemma 3.1 2 Let g a p and Q be the metric and the domain given in Theorem 3.1, and consider 
the Cauchy problem for Cf = g with /(0) defined on &o. Then, for any < s < t, there exists a 
solution in H^^iti) satisfying the energy inequality, 

11/00 HW.) < c (\\f(o)\\ 2 HttM *o) + J* \\9(T)\\ 2 H ^) dT ) > 

where C is a constant depending only on g a p, &>, and fi. 

As a result, we have two distribution functions / and f v in a common domain f2, and we can 
see that f v is non-negative by reasoning as above. For in the case of the modified equation the 
derivative of / along the particle path is strictly positive at the point (t*,a;*,j?*), a contradiction. 
Non-negativity of / is now proved by showing continuous dependence on r\ for f v . 

Non- negativity of /. We subtract the modified equation from the original one and write F = f — f v 
to get 

CF = lQ(/, /) - \Q(f n , /,) - T)e-W 2 . 

pU p u 

By applying Lemma 3.1, we obtain the following integral inequality. 

2 \ 



\\F(t)\\ 2 HliM ^<clv + l 



\Q{fJ){r)-\Q{f v J v )(r) 

pU pi) 



dr 



where we used F(0) = /(0) — f v (0) = 0. From the bilinearity of Q, the above difference of collision 
terms can be written 

Q(f, f) - Q(U fv) = 0(/. F) + Q(F, fv)- 
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We now apply Lemma 2.5 to obtain 



TOII^) <^ + /V(r)||^ T) dr) 



which implies that lim ||.F(£)||fj N (Cj t ) = 0, and consequently lim ||.F(i)|| w /fa = 0. Since [i > 5, 

we obtain 

lim||/-/ ? ,|| LM(A) =0, 

and this proves non-negativity of /. 

The above argument actually proves a stronger statement. Consider a situation where the initial 
data are not everywhere non-negative. If (t,x,p) is a point of phase space such the distribution 
function is non-negative at the point where the particle path through (t,x,p) meets the Cauchy 
surface then we can argue as before to get a contradiction. Thus a statement is obtained about the 
non-negativity of / at certain points of phase space. 

4 Further considerations 

An important concept in general relativity is that of general covariance. In the context of mathe- 
matical relativity this means that it should be possible to express the Einstein-matter equations in a 
form which is invariant under diffcomorphisms. This assertion applies in particular to the Einstein- 
Boltzmann system. The results of the previous sections are concerned with solutions of the EB 
system expressed in a local coordinate system. In other words they are results about the reduced 
Einstein equations in a harmonic coordinate system. These can be used to establish properties of 
the Einstein-Boltzmann sytem which are diffeomorphism invariant by standard procedures (cf. [16] , 
Chapter 9). From this point of view it is natural to prove a global positivity result for the EB system 
and for this it is enough to show that a solution of the Boltzmann equation on a globally hyperbolic 
spacetime with non-negative initial data is non-negative. A problem in formulating a theorem of 
this kind is to identify a suitable class of collision cross-sections. 

The conditions on the collision cross-section used in [3J and in the previous sections are conditions 
on the function S and, as such, are coordinate dependent. For a global theorem it is necessary to 
formulate an invariant condition and the function S is not invariant since it depends explicitly on 
the spacetime coordinates. A better alternative would be to formulate a condition on a. In the 
formulation of the reduced EB system using a coordinate basis the function a still depends on the 
metric components. If instead the EB system is formulated in an orthonormal frame in the way 
explained in the introduction this dependence on the metric is eliminated. This formulation also 
has an advantage for the consideration of the question of what collision cross-sections are physically 
reasonable. When the orthonormal frame approach is used the scattering kernel in general relativity 
is identical to that in special relativity. Thus the problem of identifying physically reasonable 
scattering kernels is reduced to the corresponding problem in special relativity which is better 
understood. 

What has just been said provides a strong motivation for considering the reduced system in the 
orthonormal frame formulation in more detail. The method of proof of the theorems of [3J extends 
easily to this case. The inequalities on the derivatives D% S in Lemma 12.51 are replaced by the 
analogous estimates for D*S which are actually somewhat simpler. In fact there is one estimate 
which is necessary for those proofs which is not included in 3 . In the iteration used in that paper 
there is a new metric in each step of the iteration. In estimating the differences of iterates it is 
in particular necessary to estimate the change resulting from making a change of the metric in 
the collision term. This is not mentioned in [3J. The orthonormal frame formulation eliminates 
this problem since in that case the collision term does not depend on the metric components. The 
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conditions on D^S which are analogous to those of Lemma 12.51 are coordinate independent and 
are thus appropriate for formulating a global theorem. Let us call these the orthonormal frame 
regularity conditions. From the point of view of comparing with physical conditions on the collision 
term it would be of interest to know what conditions on a are required to imply the orthonormal 
frame regularity conditions. For instance, as a simple case an assumption on the support of g can be 
made. Let us consider the orthonormal frame regularity conditions in more detail. The conditions 
of Lemma 12.51 are replaced as follows for the case N = oo. Under the same conditions on fx, k, and 
r, the following inequality should hold. 

\dia(g,e)\ 2 g 2i - i \ k '\e-^\q a ) 14 ^+ 2 ^(p°) i ^- 1 dndq < C, 

where i < \k'\ < \k\ and, motivated by the form of the expression S = Xga, the quantity a = go 
has been introduced. Suppose now that a is a smooth function and has support contained in the 
set defined by m < g < M for some positive numbers m and M. Then, the term g> 2l ~ 4 l fe I is clearly 
bounded. Moreover, the condition g < M gives the following estimate. Since we are working in an 
orthonormal frame, it implies that 

p°<\{q-p + M')<\{\q\\p\+M') for W = ]-M 2 + 1. 
q" q u 2 

After some calculations and using the inequality \p\ < p° we obtain \p\ < 2M'q a , and this implies in 
turn that p° < Cq° for some C. Hence the orthonormal frame regularity conditions hold. 

Consider a solution / of the Boltzmann equation on a globally hyperbolic spacetime and suppose 
that the scattering kernel satisifes the orthonormal frame regularity conditions. If the initial data 
are non-negative then / is non-negative. We suppose that the statement is false and obtain a 
contradiction. If the statement is false there is a point (t\,xi,pi) with f(ti,Xi,pi) < 0. Let 7 be 
the particle path passing through (ti,Xi,pi). If (t,x,p) is a point of 7 sufficiently close to the initial 
hypersurface then local coordinates can be defined such that (£, x) is contained in their domain of 
definition and the conditions of the theorem of [3] are satisfied. Hence f(t,x,p) > close to the 
initial hypersurface. Let (£2, %2,P2) be a point on 7 with /(£, x,p) > for all points on 7 with £ < £2 
and fit, x,p) < for some points on 7 with £ arbitrarily close to £2. We can choose a local Cauchy 
surface through the point (£2, £2) and a local coordinate system on a neighbourhood of (£2, X2) such 
that the conditions of the theorem of [3] are satisfied. It then follows by the result stated at the end 
of Section 3 that / > at all points on 7 with £ slightly greater than £2, a contradiction. 

It should be noted that as soon as results on well-posedness of the Cauchy problem for the 
EB system can be extended to wider classes of scattering kernels the positivity result can also be 
extended, provided the theorem includes a statement about continuous dependence of the solution 
on the scattering kernel in a suitable sense. For then it suffices to approximate a scattering kernel 
of the new class by a sequence of kernels of the class previously treated. 

We hope that by clarifying a number of issues the results of this paper will contribute to the 
development of a mature theory of the local and global Cauchy problem for the Einstein-Boltzmann 
system in the near future. 
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